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Abstract 

We make use of wavelet transform to study the multi-scale, self similar behavior and deviations thereof, 
in the stock prices of large companies, belonging to different economic sectors. The stock market returns 
exhibit multi-fractal characteristics, with some of the companies showing deviations at small and large scales. 
The fact that, the wavelets belonging to the Daubechies' (Db) basis enables one to isolate local polynomial 
trends of different degrees, plays the key role in isolating fluctuations at different scales. One of the primary 
motivations of this work is to study the emergence of the k~ 3 behavior [l[ of the fluctuations starting with 
high frequency fluctuations. We make use of Db4 and Db6 basis sets to respectively isolate local linear 
and quadratic trends at different scales in order to study the statistical characteristics of these financial 
time series. The fluctuations reveal fat tail non-Gaussian behavior, unstable periodic modulations, at finer 
scales, from which the characteristic fc~ 3 power law behavior emerges at sufficiently large scales. We further 
identify stable periodic behavior through the continuous Morlet wavelet. 

Keywords: Non-stationary Time Series, Wavelet Transform, Fractals, Power Law. 
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1. Introduction 

Study and characterization of fluctuations in financial time series are subjects of considerable interest 
0, Q (and references therein) . These non-stationary time series are well known to exhibit multiple character- 
istics depending on the scale of observation [lj. Apart from bursty behavior, fluctuations exhibit self-similar 
features, with non Gaussian characteristics [ij and inverse cubic law behavior [l], 0, [H, @|. The advent of 
wavelet transform Q enabled one to isolate fluctuations at different scales, for their characterization. The 
study of fluctuation functions, both by multifractal detrended fluctuation analysis (MFDFA) and wavelet 
analysis reveals multi- fractal behavior [5J, [6j, l8Hllj . 

In particular, the wavelets belonging to the Daubechies family enabled one to isolate local polynomial 
trends of different orders, in a manner where the fluctuations at different scales are completely independent. 
This is because of the exact orthogonality property of the discrete wavelet basis sets at different scales [l2| • 
This study has revealed non-statistical features in certain financial time series, in the context of BSE-100 
price indices The presence of correlations in various time series that can removed through shuffling is 
called the non-statistical property of a time series [llj . These non statistical correlations are the artifacts 
of external impulses provided by the socio-political factors. The fact that, determination of multi-fractal 
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behavior requires isolation of local fluctuations, both at low and high scales, motivates one to study their 
statistical characteristics more carefully. Apart from the study of their non-Gaussian properties at different 
scales, this analysis will also help in throwing more light on the non-statistical nature of the fluctuations. 
We make use of Db4 and Db6 discrete wavelets to respectively extract and study the fluctuations at multiple 
scales several large-cap equities listed in the New York stock Exchange, belonging to different sectors of the 
economy. We have systematically checked non-Gaussianity and self-similar characteristics of the multi-scale 
variations, while checking emergence of k~ 3 behavior of the averaged fluctuations. The study reveals the 
difference in the self-similar behavior of the various companies and presence of non-statistical behavior at 
small and large scales in certain cases. Fluctuation analysis indicates the presence of unstable, short periodic 
modulations. One also observes local periodic behavior of longer duration, ably identified by the Morlet 
wavelet. It is found that companies belonging to different sectors of the economy show distinct phase lags, 
in their periodic modulations. To summarize the main results of the work, in this article we show the 
convergence of the a — ¥ — 3 as a function of the scale (corresponding to coarser grained analysis) and the 
presence of periodic phase shifted modulations of the various companies as functions of different scales. We 
must point out that this kind of analysis on the fluctuations to ascertain the behavior of the power law 
exponent and the statistical properties of the fluctuations has not been carried out before to the best of our 
knowledge. 

2. Theory and Methodology 

2.1. Data 

We use the closing prices of fifteen large-cap companies listed in the New York Stock Exchange (NYSE) 
fl3j in our analysis. In this work, we limit ourselves to the behavior of the large-cap equities and are 
interested in the connections this specific section of the economy shows. The list of the equities is given in 
Table □ 



Company 


NYSE ticker 


Start date 


End Date 




symbol 


MM-DD-YYYY 


MM-DD-YYYY 


Bank of America 


BAC 


05-29-1986 


10-02-2009 


Boeing Co. 


BA 


01-02-1962 


10-02-2009 


General Electric Co. 


GE 


01-02-1962 


10-02-2009 


British Petroleum pic. 


BP 


01-03-1977 


10-02-2009 


Coca Cola Co. 


KO 


01-02-1962 


10-02-2009 


Pepsi Co. 


PEP 


01-03-1977 


10-02-2009 


Lockheed Martin Corp. 


LMT 


01-03-1977 


10-02-2009 


Walt Disney Co. 


DIS 


01-02-1962 


10-02-2009 


Wal-Mart Stores Inc. 


WMT 


08-25-1972 


10-02-2009 


International Business Machines Corp. 


IBM 


01-02-1962 


10-02-2009 


McDonalds Corp. 


MCD 


01-02-1970 


10-02-2009 


Ford Motor Corp. 


F 


01-03-1977 


10-02-2009 


Alcoa Inc. 


AA 


01-02-1962 


10-02-2009 


Merck and Co. Inc. 


MRK 


01-02-1970 


10-02-2009 


American Express Co. 


AXP 


04-01-1977 


10-02-2009 



Table 1: Summary of the data used in the analysis. The data has been obtained from |13|| . 



2.2. Wavelet transforms 

In this work, we will use both continuous and discrete wavelet transforms for the analysis of the equities. 
The Continuous Wavelet Transform (CWT) is useful in analyzing the periodic modulations present in the 
time series at different scales. It is noteworthy to point out here that since the continuous wavelet basis is 



2 



over complete, the wavelet coefficients obtained from CWT cannot be used to reconstruct the time series and 
hence cannot be used to extract fluctuations. We employ Discrete Wavelet Transforms (DWT) to achieve 
the extraction of fluctuations. 



2.2.1. Discrete Wavelet Transform (DWT) 

Discrete wavelet transform (DWT) [l2|, EH , is composed of two kernels called the father wavelet <p(t) and 
the mother wavelet ip(t), which satisfy the following admissibility conditions fh. HBj: 



(pdt < oo, / ipdt = 0, / 4>*ipdt = Q, (f) 
\0\ 2 dt = [ \^\ 2 dt = 1 , (2) 



The daughter wavelets (which are scaled and translated versions of the mother wavelet) 

ip j k (t) = Z>/ 2 ij>(2 j t -k), k <= R, j eZ+ (3) 

differ from the mother wavelet in terms of the height and the width which, at the j th scale differ by a factor 
of 2 J and 2^ 2 respectively, j and k are the scaling parameters and are real and positive integers respectively. 
The discrete wavelet decomposition of a function f(t) is given by, 

oo oo 

/(*)= E c Mt)+ J2 E d i.^i,fe(i) (4) 

k— — oo k— — ooj>0 

where c^s are the low pass coefficients, which capture the trend (average behavior or the long wavelength 
components) and dj^s are the high pass coefficients that extract the high frequency components (or fluctu- 
ations) of the function. Multi Resolution Analysis (MRA) Q leads to, 

c i. k = h ( n ~ 2fc ) c J+i<« ( 5 ) 

n 

d 3,k = ^H n - 2fc) c J + l," ( 6 ) 
n 

where, h(n) and h(n) are the low pass and high pass filter coefficients respectively. It is worth emphasizing 
that low pass and high pass coefficients at a given scale can be obtained by low pass coefficients alone at 
a higher scale. Thus, we can go from a high resolution analysis of the function to a progressively coarse 
grained picture by convolution with the filter coefficients. 

2.2.2. Continuous Wavelet Transform (CWT) 

Continuous wavelet transform was pioneered by Morlet [12j]. An integrable, well localized (in both the 
physical and Fourier domains), zero mean function, the mother wavelet ij){n), is used as the analyzing 
function. Given a discrete data set X = {x n : n £ Z + }, the wavelet coefficients are calculated by convolving 
the X, with the scaled and translated ip(n). The wavelet coefficients are given by, 

N-l / _ l\ 

Wn (s) = j2^r( !!L ^ L ) (7) 

where s is the scale. The admissibility conditions are given by, 

\m\ 2 ^ < oo, (8) 
w 

where, i>(k) = — \- f i>(x)e~^ 3 d r x (9) 
(2ir) r J Rr 

and / 4){x)d r x = 0. (10) 





(a) CWT of DIS from scale 384 to 768 



(b) CWT of KO from scale 384 to 768 



Figure 1: The scalograms representing the coefficients of the Morlet wavelets of the companies | (a) | DIS and |(b)p O, from scales 
384 to 768, the x-axis represents time, y-axis scales and z-axis wavelet coefficients W n (s). 

Periodic modulations are clearly visible. 



Here r is the number of spatial dimensions. In our analysis, we use the complex Morlet wavelet given by, 
il>(n) = Tr-VV^e-" 2 / 2 (11) 

where, n is the localized time index. Though Ea. ljlll) is a marginally admissible function, it is made 
admissible by taking loq = 6. The Fourier wavelength of ip( n )i \ is given by, 

4ns 



A 



1.03s 



(12) 



implying that the scale and the wavelength for the Morlet wavelet are approximately equal. The scale and 
frequency v are related by, 

1 



V cx 



(13) 



The cone of influence (COI) defined in jlfjj, is the region, beyond which, convolution errors make the wavelet 
coefficients unreliable for analysis 



2.3. Returns and autocorrelation function 

It is well known that for a given time series {x n }, a slowly decaying autocorrelation function given by. 

N-71 , w 

\ x n ~ Ml )\ x n+T 



S(n) 



1 



M2) 



N-n-1 ^ 
t=i 



0-1(72 



(14) 



signifies non-stationarity in {x n } [13, Here, S(n) is the correlation for the lag n € [0,37V/4], where N 
is the length {x n } and 3iV/4 is the maximum allowed lag 18]. /ii,A*2 are the mean of {x\, . . . ,xn-u} and 
{x n+ i, . . . ,xn} respectively and <J\,a-x are the corresponding standard deviations. 

However, a differencing operation, like log normalized returns can make the corresponding time series 
stationary [18j . i.e., the autocorrelation function reaches the level of noise very quickly. The normalized log 
returns (NLR), R are defined as, 

R{n) - (R(n)) 



R = 



V(R(n)*) - (R{n)) 2 



(15) 



where, R(n) — logx n +i — loga; n and J (R(n) 2 ) — (R{n)) 2 is called the volatility of returns. We consider 
the NLR for further analysis. We observe that though the stock prices show non- vanishing autocorrelation 



(Fig 2(a) I, the auto correlation function of the normalized log returns tends to zero very quickly, signifying 
stationarity(Fig 2(b) ). This behavior is observed in the time series of all the companies. We concentrate on 
the self similar characteristics of the fluctuations and use the Daubechies' family of wavelets to this end. 
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Sample Autocorrelation Function (ACF) 



Sample Autocorrelation Function (ACF) 
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(a) Autocorrelation function (ACF) of AA (b) ACF of the normalized log returns of 

AA showing zero correlation 



Figure 2: In |(a)| ones observes the non stationary behavior of a typical time series, while in |(b)| we find that the normalized 
log returns are stationary. 



2.4- Wavelet based fluctuation extraction and analysis and characterization of multifractality 
The Daubechies family of wavelets are made to satisfy vanishing moments conditions, 

t n i> jtk (t)dt = 0, (n=l,...,N) 



(16) 

due to which this family of wavelets extracts polynomial trends. For example, with Db4, where, n = 1, the 
wavelet coefficients are blind to a linear trend. These are captured by the low-pass coefficients. A discussion 
on this can be found in (tJ. We use this property to extract the fluctuations from X by reconstructing 
the original signal at each level using the approximation coefficients and then taking the difference of the 
reconstructed signal from the original signal. The algorithm used is shown in Fig.©. Fig.Q depicts the 
NLR and the fluctuations extracted by the wavelet based method. The top plot shows the NLR, the left and 
right columns are obtained using Db4 and Db6 wavelets respectively. We observe the transient and bursty 
behavior at different levels of analysis which were not apparent from the NLR. 

The discrete wavelet based method to characterize multifractality proposed by Manimaran et al. [1, Q , 
has been successfully applied to extract multifractality of various time series [l(| Ell Ell • In this procedure, 
we use the normalized log returns obtained through Ea (|15[) to obtain the time series profile by taking a 
cumulative sum of the normalized returns: 



Y(i) = ^R(t), i = l, 



,N-l. 



(17) 



Then, a DWT is applied on the profile, Y(i), to obtain the fluctuations of the profile at different scales. 
We use Db4 and Db6 to extract polynomial trends which are respectively locally linear and quadratic. The 
edge effects and the asymmetric nature affect the precision of the extracted fluctuations and hence, we use 
a DWT on the reverse profile to extract a new set of fluctuations. This new set of fluctuations are then 
reversed and averaged over the earlier obtained fluctuations to get the set of fluctuations for subsequent 
analysis. 

We subdivide the fluctuations into non-overlapping segments M s = int(iV/s), N and s being the length 
of the fluctuations and the scale respectively. The q th order fluctuation function F q (s) is then obtained by 
squaring and averaging the fluctuations over all the segments: 



F q (s) 



1 

2AL 



2M« 



(18) 
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Data: X 



NLR: R 



= R(n)-{R(n)) 
V<fl(") 2 )<fl(")> 2 



Profile: Y(i) = £* =1 R(t) 



Discrete Wavelet Transform 




Figure 3: This figure shows the algorithm for extracting fluctuations from the data. 
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Figure 4: The top plot depicts the normalized log returns of the closing stock price AA. The left column shows the fluctuations 
at levels 1, 2 and 3 respectively, using Db4 wavelet. The right column shows the corresponding results obtained by using Db6 
wavelet. 



G 



where, q is the order of the moment. We repeat the above procedure for different scales s for different values 
of q except q = 0. We thus obtain the power law scaling behavior from the fluctuation function as, 

Mi) 



F q (s) 



(19) 



in a logarithmic scale for each value of q. At q = 0, due to the divergence of the scaling exponent, a 
logarithmic averaging is employed to find the fluctuation function; 



F q (s) = cxp 



1 



2A/ 3 



2M S ^ 

3=1 



\og(F 2 (j, S )Y 



(20) 



For mono fractals time series [20j,|2l|, h is independent of q, while for multi fractals, h = h{q). The Hurst 
scaling exponent, H , a measure of the fractal nature of the time series is calculated as h(q = 2) = H and 
H 6 [0, 1], where, H g [0, 0.5) and H g (0.5, 1] reveal, respectively, the anti-persistent and persistent nature 
of the time series. 



2. 5. Tests for Gaussianity 

For a random variable x, the cumulative distribution function (CDF) is given by, 



F(x) = f 

J — oc 



f(n)dn 



(21) 



where, f(n) is the probability distribution function of X. The Kolmogorov-Smirnov test (KS-test) [22L l23j 
is used to compare F(x) of an empirical time series with that of a normal or Gaussian distribution F'{x). 
The KS-statistic is given by, 

D n =sup\F'{x)-F(x)\ (22) 

n 

where, sup is the supremum defined as the lowest clement of a subset S of a partially ordered set T which 
is greater than all other elements of S. For a continuous F(x), under the null hypothesis that it comes from 
F'(x), the quantity 



^sup\B(F(x))\ 



\fnD n 



converges to the Kolmogorov distribution given by, 
K = sup B(n) 

ne[0,l] 



(23) 



(24) 



where, B(n) is the Brownian bridge or a Brownian process. The Kolmogorov distribution is independent 
of i^(a;)[24j429|. In the present case, the KS-test indicated non-Gaussian nature for fluctuations obtained 
through wavelet transforms. 

However, we observe that though the fluctuations show non-Gaussian behavior(see Figs 5(a)|5(b)" |, suffi- 
cient averaging at the higher scales tends to bring the fluctuations towards a Gaussian behavior (skewness= 
7 = and kurtosis= k = 3), as shown in Figs 5(c) and |5(d)| 

2.6. Power law behavior of the fluctuations 

A random variable x follows a power law if it belongs to a family of distributions of the form 

p(x) cx L{x)x a , (25) 

where, L(x) is a slowly varying function such that 

L{cx) 



lim 

X^QO L(X) 



1 



(26) 



for c > 0. The It is observed (see Fig. |6(b)| that the fc -3 behavior emerges at higher scales. Also, it can 
be seen from Fig. |6(b)| that at the scale 5, all the companies seem to converge to the same value of a and 
diverge at higher scales. 



7 



(a) db4 (b) db6 (c) Skewness (7) (d) Kurtosis («) 

Figure 5: |(a)| and |(b)| depict the distribution of the fluctuations at each level against a normal distribution taken with the 
same mean and standard deviation as the fluctuations. Tl, T2 and T3 are the theoretical distributions and El, E2, E3 are the 
fluctuations at levels 1, 2 and 3 respectively. |(c)|and |(d)|show skewness (7) and kurtosis (re) against scale respectively. 



3. Multi-fractality and non-statistical behavior 

As described earlier, DbA and Db6 wavelets have been used here to extract fluctuations at multiple scales, 
from which the corresponding fluctuation functions F q (s) are constructed. The h(q) values are extracted 
from the F q (s) ~ s h ( q \ which have been plotted in FigfS] as a function of q ranging from —10 to 10. It 
is worth mentioning that the smaller values of q get dominant contributions from small functions, whereas 
the large fluctuations contribute significantly for larger q values. Multi-fractality is clearly observed in 
all the return time series. Interestingly, the /i(g)-values extracted from the shuffled time series, showed 
different behavior for some companies, particularly at low and high q- values. This reveals the non-statistical 
nature of the fluctuations of the corresponding companies j3(J ■ We carried out a Fourier spectral analysis to 
independently verify the power law behavior of the fluctuations. The cumulative of the time series revealed 
power-law behavior in the frequency domain, as is clearly seen in FiglHJ Very interestingly, one observes a 
rising component in the mid-frequency domain. This is indicative of short, unstable periodic modulations, 
as has been seen in a characteristic of chaotic dynamical system, as unstable periodic orbits. 

The dominant low-frequency periodic components can be economically extracted through the Morlct 
wavelets. The scalograms reveal two dominant modulations, depicted in FigJT] One clearly observes, phase- 
lags between the periodic modulations of companies belonging to different sectors of the economy. One also 
finds that the periodic modulation at scale 512, is non-stationary, showing varying amplitudes at different 
times. In comparison, the periodic modulations at the scale 2048, is more uniform, being present, all through 
this time domain, under investigation. The distinct phase lags physically arises from the interdependence of 
different companies with each other. This has also manifested in the study of the fluctuation characteristics 
of different companies through random matrices 31-33)]. From Fig |6(b)[ it can be observed that the k~ 



behavior emerges at higher scales. It is to be emphasised that the rising component could be related to 
unstable periodic orbits in the dynamical time series [30j ]. 



4. Conclusion 

In conclusion, we have made use of discrete and continuous wavelets to characterize the fluctuations 
in the return series of a number of large companies. As expected, the time series revealed non-Gaussian 
behavior at high-frequencies and after sufficient averaging, Gaussian behavior manifested in certain cases. 
Interestingly, the Fourier spectral analysis revealed the presence of short duration unstable modulations, 
akin to the ones observed in chaotic dynamical systems. The multi-fractal analysis showed, non-statistical 
nature of the fluctuations, for certain companies, both for high and low frequency fluctuations. It is observed 
that, after sufficient averaging the k~ 3 power law behavior emerges, when the contribution from the short, 
periodic modulations disappear. The Kolmogorov-Smirnov test performed to determine the Gaussianity of 
the fluctuations rejected the null hypothesis that the fluctuations belonged to a uniform distribution, for 
all the time scries. It was found that the skewness and kurtosis tend towards the values for a uniform 



(a) Fourier Power of IBM versus scales. 
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MRK 
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(b) Power law exponent a versus scale. 



Figure 6: |(a)| shows the Fourier power spectrum of IBM for various scales and [(b)] depicts the power law exponent (a) versus 
scale. The high frequency components reveal l//c 3 behavior. The rising component in the medium frequency range indicates 
unstable periodic modulations. 




Figure 7: The significant Morlet wavelet coefficients at scales [(a)] 512 and |(b)| 2048. The periodic behavior of the markets at 
higher scales, suggesting a long range correlation is apparent. Distinct phase lags between different companies are clearly seen. 



(a) Unshuffled returns (b) Unshuffled returns (c) Shuffled returns (d) Shuffled returns 
(Db4) (Db6) (Db4) (Db6) 



Figure 8: h(q) vs q for unshuffled and shuffled returns. 



distribution (7 s» and k = 3), as the scale increases. The Morlet wavelet identified two dominant periodic 
modulations, where different companies, showed distinct phase lags. This suggests, a systematic study of 
these periodic components, by clustering, various companies into distinct groups. Random matrix approach 
to financial time series has also yielded, distinct groups, corresponding to different sections of the economy. 
Combining wavelet based approach, with the random matrix one, will lead to a better understanding of the 
inter-relationship of different company stock values. The nature and cause of the non-statistical components 
of the fluctuations, as manifested in the difference in the h(q) values at small and large frequencies, between 
shuffled and un-shufflcd time series needs further analysis. These questions are currently under investigation. 
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